Spin waves in a skyrmion crystal 
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We derive the spectrum of low-frequency spin waves in skyrmion crystals observed recently in 
noncentrosymmetric ferromagnets. We treat the skyrmion crystal as a superposition of three helices 
whose wavevectors form an equilateral triangle. The low-frequency spin waves are Goldstone modes 
associated with displacements of skyrmions. Their dispersion is determined by the elastic properties 
of the skyrmion crystal and by the kinetic terms of the effective Lagrangian, which include both 
kinetic energy and a Berry-phase term reflecting a nontrivial topology of magnetization. The Berry- 
phase term acts like an effective magnetic field, mixing longitudinal and transverse vibrations into 
a gapped cyclotron mode and a twist wave with a quadratic dispersion. 



Baby skyrmions are magnetic textures conjectured to exist in two-dimensional Heisenberg ferromagnetsii Like 
domain walls and vortices, skyrmions are stable for topological reasons. A localized magnetic texture, parametrized 
by the unit vector field m(r), has a quantized topological charge 



A state with a single skyrmion (n = ±1) cannot be continuously deformed into a uniform ground state (n = 0). While 
skyrmions are metastable excitations in the pure Heisenberg model, the presence of additional interactions may lower 
their energy cost and create a ground state with a finite concentration of skyrmions. Bogdanov and collaborators 
pointed out that a Lifshitz invariant m • (V x rh) in the free energy may stabilize skyrmion-like line defects in a 
three-dimensional ferromagnet with easy-axis anisotropyi^^— Such a term is allowed in a crystal without the inversion 
symmetry. 

Recent experiments^^— have revealed ordered arrays of skyrmion lines in three-dimensional MnSi and Fei_a;Coa;Si 
and arrays of skyrmions in thin films of Feo.5Coo.5Si and FeGe. A common theme of these materials is the lack 
of the inversion symmetry. Thin films are a particularly favorable environment for "skyrmion crystals" (SkX) . The 
stabilizing factor is likely dipolar interactions, which tend to favor inhomogeneous states such as the closely related 
bubble phase in a two-dimensional Ising ferromagneti^ 

The dynamics of magnetization in the SkX phase is an interesting question in its own right. Zang et al. approached 
it by treating skyrmions as point-like particlesiiS Because of a nontrivial topology of a skyrmion, its kinematics are 
similar to that of a charged particle in a uniform magnetic field. Moving skyrmions are subject to a Lorentz-like force 
that is proportional to the skyrmion velocity and is thus more important than inertia in the limit of slow motion. (A 



I. INTRODUCTION 
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FIG. 1: Distribution of magnetization M(r) of a skyrmion crystal in a plane z — const with magnetic field applied normal 
to that plane. The crystal is a superposition of three helices with wavenumbers given by Eq. ([2]). The vector plot shows 
components Mx and My; shades of gray encode Mz (arbitrary units). 
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similar effect is well known for magnetic vortices?ii) Like a Wigner crystal in magnetic field )^^'^'^ an SkX exhibits a 
quadratic phonon dispersion, to ~ /2m. 

In the SkX phase, skyrmions form a dense lattice with the distance between skyrmions comparable to their size.— >^ 
In this case, treating skyrmions as point particles may not be justified. From a complementary perspective, the SkX 
phase is a spin-density wave (SDW) depicted in Fig. [T] It is a superposition of uniform magnetization parallel to the 
applied field H and three phase-locked helices whose wavevectors, 

koa = A:o(l,0,0), ko6 = fco(-l/2,\/3/2,0), koc = fco(-l/2, -%/3/2, 0). (2) 

form an equilateral triangle in the plane normal to H.^ In the vicinity of the Curie point, higher harmonics are 
suppressed because they are not soft modes. This makes the SDW picture a more accurate starting point. In light of 
that, we set out to characterize low-frequency excitations of the SkX phase in the SDW limit. The calculation of the 
spin-wave spectrum is nontrivial even for a sing le helix.iiii^ To make the paper self-contained, we present a simple 
derivation of the spin- wave spectrum in a single helix for a three-dimensional ferromagnet in Sec. In] We then discuss 
the case of three coupled helices in Sec. IIIIl 

The results can be briefly summarized as follows. The SkX (Fig. [T]) breaks the symmetry of translations in the 
xy plane; translations in the z direction remain a good symmetry. Thus low-energy excitations of the crystal are 
associated with displacements in the xy plane, u(r) = {uxTUy,Q). The potential energy of these displacements is 
expected to be the same as that of an isotropic two-dimensional (columnar) solid, 

V- V- /^^ \ B{dlnY 

2^ 2^ {-u^iUjj + ^lUijUij \ + . (3) 

i—x^yj—x,y ^ ^ 

Here Uij = {diUj + djUi)/2 are components of the strain tensor in the xy plane, and A and /i are the elastic Lame 
coefficients.^^ Displacements inhomogeneous along the z-axis bend the crystal, hence an elastic energy proportional 
to (d'^u)^ with a bending modulus B. The dispersion of elastic waves depends on the kinetic terms in the Lagrangian. 
By analogy with ordinary solids, one might expect a kinetic-energy term p(u)^/2, which is indeed present. Excitations 
in a columnar solid with these properties would be longitudinal and transverse sound waves with a linear dispersion 
for waves propagating in the plane of the crystal and a quadratic one for waves propagating normal to the plane. 
However, the nontrivial topology of magnetization in a unit cell adds a Berry-phase term Bu^Uy to the Lagrangian of 
the SkX. Because this term is linear in the velocity, it dominates over the kinetic energy in the limit of slow motion. 
The resulting Lagrangian is 

r r,- V- V- \ B{dlvLf 

C = Bu^Uy + — 2^ 2^ [l^UuUjj + ^lUiJU^■ij , (4) 

i—x^yj—x^y 

At the lowest frequencies, transverse and longitudinal spin waves are mixed and exhibit a quadratic dispersion, 
Lo ^ /2m, in agreement with the result obtained in the limit of point- like skyrmionsiiS 

With the sole exception of the Berry-phase term, the coupling constants of the SkX are set by the physics of a 
single helix. It is interesting that the anharmonic coupling that locks the phases of the helices does not enter the 
Lagrangian Q. This coupling constant sets an upper limit of frequencies for which the three locked helices move 
together as a solid. 



II. SINGLE HELIX 



A. Static solutions 



We consider an isotropic ferromagnet in three dimensions with the magnetization field M(r). The Landau free- 
energy functional is 

F[M(r)] = j (fr [^(VM)^ + aM^ + cM'^ + DM ■ (V x M)] , (5) 

where A is the exchange coupling and (VM)'^ = J^i ■ The Dzyaloshinskii-Moriya coupling DM ■ (V x M) 

breaks the inversion symmetry and is only allowed in magnets with a noncentrosymmetric lattice. 
The quadratic part of the free energy is diagonalized by Fourier transform, 

F[M{r)]=Y,MQ*AfMl, (6) 

k 
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FIG. 2: The global frame {fii, 112,113} and the local frame {61,62,63} of a single helix, Eq. 



with a coupling matrix 



A^'' = (a + - iDe^'^Tp. 



(7) 



Its lowest eigenvalue, a + Ak^ — Dk, is minimized by wavenumber fco = D/2A. When a drops below the critical value 
flc ~ D"^ /AA, the ground state is a helical SDW, 



M(r) = Mo[ni cos (ko • r) + 112 sin (ko • r)], 



(8) 



where ko — koh^ and A/g = (ac — a) /{2c). The unit vectors {£11,112,113} are mutually orthogonal and form a 
right-handed triplet. We will also find it convenient to use a local frame set by the unit vectors {61,62,63}, where 



ei = 111 cos (ko • r) + n2 sin (ko • r), 62 = -fii sin (kg • r) + n2 cos (ko • r), 63 = 113. 



(9) 



Here ei is parallel to the local direction of magnetization in equilibrium, the pair (61,62) defines the plane of the 
helix, and 63 is parallel to its wavevector ko. The orientations of the local and global frames are shown in Fig. [2] 



B. Spin waves 



Small deviations from equilibrium can be described by two angles parametrizing twists of magnetization in the 
plane of the helix (a) and out of the plane {/3): 



M/Mo = 61 [1 - {a^ + p^)/2] + ae2 + ^63. 



(10) 



For convenience, we align the coordinate axes x, y, and z with the vectors fii, n2, and 113. To the second order in the 
twist angles, the energy cost is 



U = DMqI—Ox sinkoz + dy cos kQz)l3 
+AM^ [(Va)2 + (V/3)2 + kll3^] 

+DAfo [/^(Sa; COS koz + dy sin koz)a — a{dx cos koz + dy sin koz) 



(11a) 
(lib) 
(11c) 



The linear term (jllap is a total divergence that reduces to a boundary term upon integration over the volume. It 
does not influence the physics in the bulk and can be ignored. 

The dynamics of spin waves can be obtained by constructing a Lagrangian for the fields a and /?. The kinetic 
term comes from the Berry phase for a spin, Sb = S J (cos 6* — l)(pdt.^''^ Measuring the inclination angle 9 from the 
local magnetization axis ei yields a = 9cos(f) and /? — 9sin(f) for infinitesimal a and /3, so that the Berry-phase term 
in the Lagrangian becomes Cb — J{a.P — a/3)/2, or equivalently Ja.p, where J = Afo/7 is the density of angular 
momentum and 7 is the gyromagnetic ratio. The resulting Lagrangian is 



L = Lb-U = JaP - AMI [i^oi? + f^f + ^0/^^] - 2DM^l3{dx cos koz + dy sin kQz)a. 



(12) 
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The Euler-Lagrange equations of motion are 

- [Jdt + 2DM^{d^ cos koz + dy sin koz)](3 = 2AM^\/^a, (13a) 
[Jdt + 2DM^{d^coskoz + dySmkoz)]a = 2AM^{V^ P ~ k^^). (13b) 

The equations, though linear, are comphcated by the presence of oscihatory terms reflecting the breaking of the 
translational symmetry by the helix. Simple closed-form solutions can only be obtained for waves propagating along 
the direction of the helix h^, a = a cos {ujt — q^^z), P — jS sin {ujt — q\\z), with the dispersion 

= i2AMS/J)^k^,q^^+qf^. (14) 

For short wavelengths, q\\ 3> ko, we recover the quadratic dispersion of magnons in a Heisenberg ferromagnet. For 
wavelengths longer than the pitch of the helix, the dispersion becomes linear, uj ^ sg||, with the wave velocity 
s = 2AM§kQ/J = jMqD. It is worth noting that the wave velocity is independent of the exchange coupling A. 



C. Long- wavelength approximation 

Further progress can be made by analyzing spin waves in the long- wavelength limit. The helical modulation mixes 
waves with wavevectors q and q ± ko, but the effects of the mixing can be controlled in the limit g <C fcg. 

As usual, we expect that low-frequency waves are associated with the generators of broken symmetries, in this case 
the global rotational symmetry SO (3) and translations along the helix axis 113. Zero modes, which leave the energy 
of the system invariant, can also be obtained by requiring that the first variation of energy vanish: S J U dV = 0. The 
general solution is an infinitesimal global rotation parametrized by three angles {(t>i} about the global axes {Ai}: 

a — ~(l)ikQy + (l)2kQX + (1)3, /S — (pi sin koz ~ (1)2 cos koz. (15) 

Note that 03 can also be viewed as a translation of the helix through distance —(fis/ko along its axis £13. 

Although there are three parameters describing an infinitesimal global rotation, only one of them — (j)^ — is associated 
with a Goldstone mode. Radzihovsky and Lubensky^^ pointed out that the other two modes become "massive" via 
a Higgs-like mechanism similar to that in smectics A.— 

The explicit form of the zero modes (|15l) shows that soft modes of the helix are situated near wavevectors and 
ko. We therefore focus on the respective Fourier components of a and /3: 

a(r) = ao -1- ai cos fcoz + a2 sin fcoz, l3{r) — /Sq + /3 1 cos k^z + (32 sin kg z, (16) 

where fields ai(r) and /3i(r) vary slowly in space. We express the Lagrangian (fT2|) in terms of these. After averaging 
over spatial regions large compared to the helix period, oscillatory terms such as cos2kQZ vanish and we obtain 

£ = do/3o - (Vao + /3)2 - /?2 _ IZ^I! (17a) 

+ ^ - (V/3o)' - ^ - + 2a • V/3o + ^ + ^ ■ (17b) 

Here we introduced vectors a. — (ai, a2, 0) and (3 = (/3i, /32, 0) and switched to natural units, 

length: 1/fco, time: J/AMlkl, energy: AMg/fco. (18) 

The term {V jSf is understood as {V Pif + (V/32)^. 

The low- frequency dynamics of the soft mode ao is determined by the terms displayed in Eq. (|17al) . The rest 
of the terms (I17bp provide corrections containing higher powers of the gradients and can be safely dropped in the 
long- wavelength limit. 

Integrating out the hard field /3o with the aid of its equation of motion, /3o — <io/2, generates a "kinetic energy" for 
the soft mode, q;o/4. At low frequencies and long wavelengths, the "gauge field" (3 follows the gradient of the "Higgs 
field" ao so that /3 w — V^ao, where Vj_ = {dx,dy, 0) is the transverse part of the gradient. Integrating out (3 yields 
the following Lagrangian for the soft mode ao: 

/: = ^-(v,|ao)^-^. (19) 
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FIG. 3: Spin waves in a helical state. Left: a helix in equilibrium. The helix wavevector ko is vertical. Solid red and dashed 
green lines are planes of constant phase koz = 27rn and koz = 27r(n + 1/2). Center: a spin wave with the propagation vector q 
parallel to ko. Right: a spin wave with the propagation vector q perpendicular to ko. Planes of constant phase have a vertical 
displacement given by Eq. (|20|l . 



Here V|| = dz is the longitudinal part of the gradient. 

Low- frequency spin waves are slow twists of the magnetization in the helix plane described by the "Higgs field" ao, 
Fig. [3] Twists out of the plane, described by the "gauge field" /3, adjust to them as needed. Alternatively, the spin 
waves can be viewed as displacements of planes of constant phase of the helix along its propagation direction 63, 

u{r) = -ao(r)/fco. (20) 

Waves propagating along the helix direction are regular phonons, with the energy density (V||ao)^- Waves propagating 
along the planes of constant phase are fiexural modes, with a much lower energy density (V^ao)^/2. This results in 
a strongly anisotropic dispersion. With the units of length and time restored, we have 

u; = {2AMllJ)^klql+qy2. (21) 

This expression agrees with the dispersion for the longitudinal direction, Eq. ([Til) , in the long- wavelength limit. 



III. SKYRMION CRYSTAL 



A. Three coupled helices 



The skyrmion-crystal phase is realized in the vicinity of the critical point in an applied magnetic field H.^ The field 
induces a finite uniform magnetization M — xH. In its presence, the quartic term in the free energy ([5]) generates an 
effective cubic term, 

= 4c(M-M)(M-M), (22) 

where the order parameter M contains only soft modes with wavenumbers near fco and excludes the uniform component 
of magnetization. By the standard argument (^li^i^S the cubic term favors a magnetic state that is a superposition of 
three helices ^ whose wavevectors koa, kob, and koc form an equilateral triangle orthogonal to the applied field H. 
(It should be noted that the helix amplitude Mq need not be the same as in the case of a single helix. Here we treat 
Mq as a phenomenological parameter that determines the energetics of the soft modes.) The cubic term couples the 
phases of the three helices thereby generating an eff'ective coupling between their soft modes aoa, ao6, and a^c'- 

II / , , \ , cr(Q!Oa + "Ob + aoc)^ .„„^ 

"lock = -cr COS {oLQa + aob + "Oc) ~ H . (23 ) 



where 



_ %cMMl _ 4M 

- ~AMi^ - iWofcleF ^^^^ 
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is a dimensionless phase-locking coupling and = \J Aj (0,^ — a) is the longitudinal correlation length of 
magnetization. ^ ^ 

At the harmonic level, the three soft modes have the following Lagrangian: 

^ ^ _aiao. + ao, + ao^r + ^ (^^ _ (V,,, a,,f iZl^^ , (25) 

i—a^h.c 

where V||j = ko; • V and V^; = kg^ x V. Thanks to the phase-locking coupling ct, one of the modes acquires a finite 
frequency = \/6o'. The two remaining Goldstone modes are translations u in the plane of the SkX. (Translations 
normal to the plane do not alter the state of the system.) For convenience, we direct the z-axis along the applied field 
H and choose the x and y-axes in such a way that 

koa = (1,0,0), kob = (-1/2,73/2,0), koc = (-l/2,-\/3/2,0). (26) 

As in the case of a single helix, Eq. ([20)) . a translation u yields a phase shift a^i = — koj • u for helix i = a,b, c. The 
resulting Lagrangian for the Goldstone modes and Uy describes a columnar solid: 

P(u)2 V- V- , \ Bid^uY 



i=x,y ]=x,y 



C = — 2^ 2^ [-u,iU„+^iu,JUij ] ^ , (27) 



where Uij = {diUj + djUi)/2 are components of strain in the xy plane, p — 3/4 is the density, B = 3/2 is the bending 
modulus, and A = /i = 3/4 are the Lame parameters. Note that the phase- locking coupling a does not enter the 
Lagrangian ((27l) . Instead, it sets the range of frequencies, 

UJ < V6a, (28) 

in which the texture behaves as a two-dimensional solid. At higher frequencies, the three helices decouple from one 
another. 

Eigenmodes of the Lagrangian (|27p come in two flavors: longitudinal waves, in which the displacement of the SkX 
u ~ {ux, Uy, 0) is parallel to the in-plane component of the wavevector qin = {qx, Qy, 0), and transverse waves with u 
along z X qin = (—91;, Qx, 0). With the physical units restored, the dispersions are 

uJi = {AM^/J)^3P,qf^ + 2qt^„ u;t ^ {AM§ / J) ^ k^ql + 2qt^„ (29) 

Here gout = Qz is the out-of-plane component of the wavevector. This dispersion was obtained recently by Kirkpatrick, 
Belitz, and collaborators, ^^'■^■^ who treated the SkX as a two-dimensional solid (|27| . 

The first term in the effective Lagrangian ()27|) is proportional to the square of the velocity u and thus can be viewed 
as kinetic energy of the SkX. As a result of that, waves propagating in the plane of the crystal have a linear dispersion 
u! — sq with the speed si — AM^koy/S/J' — for longitudinal waves and sj ~ s//\/3 for transverse ones. 



B. Berry phase 

The low- frequency dispersions will be modified if the effective Lagrangian (|27p acquires a Berry-phase term BiixUy. 
Being linear in the velocity, this term dominates over kinetic energy in the limit w — > 0. The effective Lagrangian 
then becomes 

. „. p(u)' V- V- /^A \ B(a2u)2 

C = BuxUy + — 2^ 2^ (2"""^^ I ^ • (30) 

i—x.yj—x^y ^ ' 

We focus on waves propagating in the plane of the crystal, so that S^u = 0. The Berry-phase term acts like an effective 
magnetic field hybridizing longitudinal and transverse waves. One of the resulting modes has a finite frequency, 

= B/p, (31) 

in the g — >■ limit. It is the cyclotron mode that reflects the interplay of inertia and the Lorentz force acting on a 
charged particle in a magnetic field. The other mode has a quadratic dispersion, 

w = ujiutIB = {qlJB)^Jp{\ + 2p). (32) 
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This result holds for sufficiently small in-plane wavenumbers, gin ^ B. In the opposite limit, gin ^ B, the motion 
becomes sufficiently fast for the Berry phase to be negligible relative to kinetic energy and the eigenmodes become 
longitudinal and transverse sound waves. 

The Berry-phase coupling B can be computed by moving the SkX along an infinitesimal contour of area AX AY: 

u u + (AX, 0,0) + {AX, AY, 0)^u+ (0, AY, 0) ^ u. (33) 

In the process, the magnetization vector at a given point traces out a closed path on a sphere: 

dM dM dM dM 

M^M-^AX ^M-^AX -^AY ^M-^AY ^M. (34) 
ox ox oy Oy 

If magnetization is not confined to a single plane then the path on the sphere has a nonzero area and the spins acquire 
a Berry phase Sb proportional to the area on the sphere and spin length. In this way we obtain 



„ I dx dv M ^ , ^ ^ ^ ^ ^ , ^ 

B = - I m • (a,m X dyxh) , (35) 



where rh = M/M is the unit vector parallel to magnetization and the integration is done over a unit cell of area 

Deep in the ordered phase, the magnetization length M is fixed and Eq. ([55]) reduces to 

AnnM , , 

^ = — ^' 

where 

^^^^m^^ni^^ (37) 
47r 

is a topological charge of the unit cell known as the skyrmion number. In an SkX, n = ±1. In natural units, 

B^^^. (38) 
27rMoV3 

Closer to the critical point, an SkX is a superposition of three phase- locked helices on top of uniform magnetization. 
In such a texture, magnetization length varies in spacei^ Then nM in Eqs. ([55)) and ([55]) should be understood as 
the average length of magnetization weighted with skyrmion density rh • (9a;m x S^rii) /(47r). 



IV. DISCUSSION 



We have derived a low-energy description of a skyrmion crystal viewed as a superposition of three phase-locked 
helices.^ This approach is justified because skyrmion crystals are typically observed near the critical point, where 
soft components of magnetization have characteristic wavenumbers fco- In this limit, the extent of a skyrmion is 
comparable to the interskyrmion separation, which means that these solitons cannot be treated as point particles. 
Viewing this magnetic texture as a superposition of three helical waves is a better starting point. 

Our main result is the Lagrangian for low-energy excitations of the skyrmion crystal, Eq. pop . It is expressed in 
terms of the Goldstone modes of a columnar solid in three dimensions, deformations parametrized by a two-component 
displacement field u(r) = (ita:,Wy,0). The description is valid at sufficiently low frequencies, when the three helices 
behave as a cohesive solid. The upper limit is set by the strength of the anharmonic coupling constant locking 
the phases of the helices. At intermediate frequencies, the dynamics is dominated by the inertia of the crystal, 
parametrized by the kinetic energy in Eq. (1301) . In this regime, excitations are transverse and longitudinal vibrations 
with a linear dispersion, cj ^ sq, for waves propagating in the plane of the crystal. 

At the lowest frequencies, the Berry-phase coupling mixes longitudinal and transverse waves. One of the mixed 
modes acquires a quadratic dispersion, a; ~ g^/2m p2p. as proposed by Zang et al}^- The other mode has a ffiiite 
frequency (j3ip in the g — >■ limit. This mode corresponds to cyclotron motion of skyrmions: the Lagrangian (j30p 
indicates that skyrmions behave as massive particles with density p in an effective magnetic field of strength B. The 
skyrmion mass density in standard units is p = 3J'^/(4^Mq). The mass of a skyrmion nis can be compared to the 
mass of a magnon mm in the same ferromagnetic material (without the Dzyaloshinskii-Moriya coupling): 



mm a {koay ' 
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where £ is the length of the skyrmion Une (in a three-dimensional crystal) and a is the lattice spacing of the material. 

The observability of the skyrmion cyclotron mode depends on the strength of the skyrmion crystal measured by the 
phase- locking coupling a (|24|) . The cyclotron should be in the range of frequencies (f28| . where the texture behaves as 
a two-dimensional crystal, rather than three independent helices. If the crystal is too weak, a ^ the Lagrangian 

(pO| only applies in the limit of slow motion, where inertia of the skyrmion crystal can be neglected. In that case, the 
density term may be entirely omitted. 

It is interesting to compare the spectra of low-energy excitations of skyrmion crystals in ferromagnets and in systems 
exhibiting the quantum Hall effect (QHE).^^ In QHE systems, the Coulomb interaction between charged skyrmions 
in two dimensions causes the Lame parameter A to diverge as at low wavenumbers, which leads to the g'^/^ 
magnetophonon dispersioui^i^i The cyclotron mode ([5T|) is absent because the mass density p is smalli^ This is 
similar to our case when a small p pushes the cyclotron frequency outside the frequency range (pS)) in which the 
texture can be thought of as a skyrmion crystal. In the absence of spin-orbital coupling, QHE skyrmion crystals get 
an additional Goldstone mode with a linear dispersion, associated with the spontaneous breaking of the global SO (2) 
symmetry of spin rotations.— In our system, such spin waves are coupled to translations of the skyrmions and are 
therefore gapped via the Higgs mechanism.^* 

Although in this paper we have focused on a three-dimensional ferromagnet,'^'^ these results can be readily extended 
to thin films,— 1^ in which long-range dipolar interactions make the SkX phase particularly stable. 

Note added in proof 

A recent numerical study by Mochizuki2£ provides evidence for two branches of spin waves in a skyrmion crystal. 
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